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ABSTRACT.  Simulations  were  made  of  the  recovery 
of  mean  gravity  anomalies  from  intersatellite 
Doppler  measurements.  A  pair  of  surface-force 
compensated  satellites  was  assumed  to  be  in 
identical  polar  orbits,  spaced  by  3°,  at  an  alti- 
tude of  200  km,  taking  intersatellite  data  with  a 
precision  of  0.001  mm/s  .  The  results  for  a 
half-year  mission  indicate  that  l°xl°  mean  gravity 
anomalies  can  be  obtained  to  a  precision  of  a  few 
milligals.  A  so-called  high-low  case,  when  one  of 
the  satellites  is  geosynchronous,  appears  to  per- 
form about  twice  as  well,  but  is  probably  imprac- 
tical to  use  because  of  the  difficulty  in  obtain- 
ing the  required  data  precision  (10~3  mm/s  )  with 
the  long  (~36,000  km)  intersatellite  distance 
involved.  In  contrast  to  previous  analyses,  our 
results  are  indicative  of  data  obtainable  from  a 
global  solution  free  of  debatable  a  priori  con- 
straints or  assumptions  about  the  geopotential . 
Thus  a  properly  optimized  mission  using  additional 
information,  such  as  results  from  altimetry  or 
mean  surface  anomalies,  could  improve  upon  the 
results  presented  here. 

INTRODUCTION 


The  concept  of  using  intersatellite  measurements  to  determine  the 
geopotential  is  at  least  10  years  old  (Wolff  1969).  The  advantage  of  this 
scheme,  compared  to  those  based  on  ground-based  observations,  is  continuous 
global  observation  of  orbital  perturbations  produced  by  the  anomalous 
geopotential. 


The  long  wavelength  variations  of  the  geopotential,  e.g.,  n  <  6,  have  been 
determined  satisfactorily  from  conventional  satellite  tracking  data.  Pro- 
gress also  has  been  made  at  very  short  wavelengths  using  GEOS-3  altimeter 
data.  Rapp  (1979)  reported  l°xl°  mean  gravity  anomalies  from  these  data 
precise  to  ~7  mgal  over  most  of  the  ocean  surface.  However,  a  global 
determination  of  the  gravity  field  to  a  higher  precision  would  seem  to 
depend  on  a  new  satellite  system. 

For  a  believable  simulation  of  the  recovery  of  gravity  field  parameters, 
the  following  conditions  are  necessary: 

•  The  simulation  should  have  a  global  character. 

•  The  simulation  should  consider  the  effect  of  random  and  systematic 
errors . 

•  The  simulation  should  be  independent  of  a  priori  information. 

The  last  requirement  may  appear  stringent,  but  adherence  produces  conserva- 
tive conclusions. 

REPRESENTATION  OF  THE  EARTH'S  GRAVITY  FIELD 

Spherical  harmonics  were  unsuitable  for  our  analyses  because  a  complete 
set  of  coefficients  would  have  been  required  to  model  purely  local  effects. 
A  spherical  harmonic  model  of  degree  and  order  180  (l°xl°  resolution)  would 
require  32,400  coefficients.  In  addition,  the  computation  of  harmonics  above 
degree  and  order  70  requires  special  care  (Chovitz  et  al.  1973).  We  chose 
instead  to  represent  l°xl°  mean  gravity  anomalies  by  a  surface  layer  of  vari- 
able density.  In  particular,  we  assumed  the  density  X  to  t>e  constant  over  a 
small  block,  and  represented  the  total  field  in  terms  of  an  aggregation  of 
such  blocks. 

Another  representation,  spherical  harmonic  sampling  functions,  has  been 
proposed  by  Giacaglia  and  Lundquist  (1972).  However,  these  have  severe 
computational  problems.   Sampling  functions  for  a  line,  plane,  or  any 


Cartesian  Space  R,  can  be  generated  by  translation;  i.e.,  one  needs  only  to 
shift  the  argument  by  one  division.   For  example,  if  we  define 

f   .  _  sin  (27tx) 
Vx;  ~    27IX 

we  can  derive  the  remainder  of  the  sampling  functions  from 

SR  (x)  =  Sq  (x  -  k),   k  ±1,  ±2,  .  .  .  . 

Unfortunately,  a  simple  symmetry  condition  does  not  exist  for  the  spherical 
domain,  so  one  would  have  to  invert  32,400  nonsparse  matrices  of  dimension 
32,400  x  32,400  to  obtain  the  functions,  which  would  be  represented  by  a 
series  with  32,400  terms.  Thus,  technical  and  economic  considerations 
determined  the  use  of  a  surface  density  model  in  these  simulations.  In  ad- 
dition, the  ORAN  (orbit  analysis)  computer  program,  described  later,  did  not 
have  to  be  modified  and  retested  extensively  before  using  the  density  model 
subroutines.  The  logic  in  Morrison's  (1977)  program  was  changed  to  accept 
blocks  of  arbitrary  size  (Morrison  1979).  We  did  not  use  more  sophisticated 
means  of  analyzing  the  results,  such  as  inversion  theory  (Kaula  et  al.  1978) 
or  collocation  (Moritz  1978)  because  the  project  was  designed  to  determine 
the  sensitivity  and  resolution  unique  to  the  proposed  satellite-to-satellite 
tracking  systems. 

The  density  of  a  block  is  related  to  the  usual  mean  gravity  anomaly  Ag 
(Orlin  1959) 

G*  =  I  +  I  I  N  -  m  +  °-0368  N 

where  G  is  the  gravitational  constant,  y  is  global  mean  gravity,  R  is  mean 
Earth  radius,  N  is  geoid  height  (m) ,  x  is  surface  density,  and  Ag  and  Gx  are 
in  milligals  .   For  most  of  the  Earth,  the  approximation 

Ag  =  2n  G  x 

suffices  for  the  purpose  of  error  analysis. 


Computing  the  acceleration  of  a  satellite  resulting  from  a  l°xl°  density 
layer  block  involves  evaluation  of  a  double  integral.  Some  authors 
approximate  the  solution  of  this  integral  by  the  attraction  of  one  or  more 
point  masses  (Koch  and  Witte  1971).  We  used  the  solution  of  Morrison  (1976, 
1977).  He  evaluates  one  integral  approximately  analytically  and  the  other 
integral  numerically.  The  advantage  of  this  latter  scheme  is  that  the  block 
has  no  artificial  "structure."  Unless  a  very  large  number  of  point  masses 
is  used  to  approximate  an  individual  density  block,  this  structure  creates 
rapidly  oscillating  perturbations  on  a  low  altitude  orbit,  possibly  leading 
to  erroneous  conclusions. 

The  global  rms  (root  mean  square)  mean  gravity  anomaly  in  1°  squares  is 
about  30  mgal.  Because  the  uncertainty  of  l°xl°  mean  anomalies  visible  to 
GE0S-3  over  the  ocean  surface  is  less  than  10  mgal,  let  us  consider  the 
effect  of  a  10-mgal,  l°xl°  block  on  the  speed  of  a  near-Earth  satellite. 
Figure  1  shows  the  velocity  perturbation  in  the  along-track  direction  for  a 
satellite  in  a  circular  polar  orbit  passing  directly  over  the  block  at  an 
altitude  of  200  km. 

Several  interesting  features  appear  in  figure  1.  First,  note  that  the 
block  has  two  effects,  a  small  high  frequency  perturbation  of  about 
0.006  cm/s,  and  a  much  larger  low  frequency  (once/revolution)  effect.  This 
latter  effect  comes  from  the  overall  perturbation  of  the  energy,  and  hence 
period  (Hotine  and  Morrison  1969).  The  effect  does  not  contribute  informa- 
tion because  every  block  will  contain  the  effect  and  adjacent  blocks  cannot 
be  separable  on  this  basis  since  the  wavelength  of  the  effect  is  vastly 
greater  than  the  dimensions  of  the  blocks.  The  significant  information  about 
a  block  appears  in  the  sharp  high  frequency  oscillation  which  occurs  as  the 
satellite  passes  over  a  block.  This  oscillation  persists,  however,  for  a 
considerable  time  (2-3  minutes)  compared  to  the  time  required  for  a  satel- 
lite to  traverse  the  l°xl°  block  (0.25  minute).  The  fact  that  the  perturba- 
tion of  a  block  persists  over  an  along-track  distance  of  about  10  times  the 
block  size  also  means  that  adjacent  l°xl°  blocks  will  have  a  similar  effect. 
These  will  be  hard  to  separate  even  when  using  the  high  frequency  impulse. 
It  is  also  apparent  that  to  maximize  the  relative  along-track  velocity  per- 
turbation between  two  satellites  in  identical  low  orbits  the  separation  of 
the  satellites  must  be  greater  than  1°. 
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Figure  l.--The  perturbation  in  the  along-track  component  of  velocity  result- 
ing from  a  10-mgal,  l°xl°  mean  anomaly  for  a  satellite  at  200-km  altitude, 


The  magnitude  of  the  velocity  perturbation  could  be  increased  by  choosing 
a  lower  altitude,  but  200  km  is  a  reasonable  height  for  a  practical  mission 
after  considering  atmospheric  drag  forces.  At  this  altitude  the  atmospheric 
density  scale  height  is  only  about  35  km,  so  the  drag  force  is  very  sensi- 
tive to  small  changes  in  altitude.  In  the  discussion  that  follows,  we 
assume  all  satellites  used  surface-force  compensation  devices  to  eliminate 
drag  and  radiation  pressure  effects. 

Figure  2  shows  the  relative  speed  (i.e.,  range-rate)  between  two  satel- 
lites passing  over  the  10-mgal,  l°xl°  anomaly  (used  in  fig.  1)  for  various 
separations  of  the  satellites.  Larger  separations  yield  a  larger  signal. 
Figure  3  shows  the  maximum  signal  occurs  for  a  separation  exceeding  5°.  It 
is  important  to  remember  that  the  maximum  signal  is  dependent  on  the 
altitude  as  well  as  the  separation,  causing  figure  3  to  be  valid  only  for  a 
200-km  altitude. 

Figures  2  and  3  illustrate  the  essential  problem  of  determining  the  fine 
structure  of  the  gravity  field  from  satellite  perturbations.  The  effect  of 
a  l°xl°  block  is  very  small.  Thus,  attempting  to  estimate  the  geopotential 
in  terms  of  l°xl°  anomalies  requires  highly  accurate  data  and  attention  to 
the  effect  of  systematic  errors. 

SIMULATION  METHOD 

The  most  difficult  aspect  in  simulating  the  recovery  of  the  global  geopo- 
tential from  satellite-to-satellite  tracking  is  the  need  for  a  global 
analysis.  In  our  situation,  we  are  interested  in  a  resolution  of  l°xl°,  but 
there  are  41,000  l°xl°  equal-area  blocks  on  the  sphere.  Solutions  proposed 
by  other  authors  have  not  been  entirely  satisfactory.  Argientero  et  al. 
(1978)  simulated  the  adjustment  of  a  small  region  of  l°xl°  blocks  and 
carried  the  additional  surrounding  blocks  as  unadjusted  parameters  with  a 
given  variance  in  the  solution.  Schwarz  (1971)  also  computed  the  covariance 
matrix  of  adjusted  blocks  in  a  limited  region,  but  with  the  tacit  assumption 
that  the  blocks  which  were  not  adjusted  were  known  without  error.  Clearly 
the  former  approach  leads  to  pessimistic  results,  the  latter  to  optimistic 
conclusions. 
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Figure  2. — Intersatellite  relative  speed  (p)  for  various  separations 

of  the  satellites. 
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Figure  3. --Maximum  variation  of  intersatellite  range-rate  (p)  for 
various  satellite  separations  (Au) . 


Our  approach  for  simulating  a  global  solution  was  to  simulate  the  recovery 
of  a  region  with  dimensions  large  enough  to  allow  the  l°xl°  anomalies  near 
the  center  of  the  region  to  be  regarded  as  having  essentially  the  same  sta- 
tistical properties  in  the  solution  that  they  would  have  possessed  in  a 
global  solution.  Figure  1  shows  the  size  of  this  region  would  have  to  be 
substantial  because  the  impulse-like  perturbation  of  a  l°xl°  block  persists 
for  a  long  time  relative  to  the  size  of  the  block.  We  began  our  experiments 
with  a  15°xl5°  region  containing  225  l°xl°  mean  gravity  anomalies.  To 
verify  that  this  region  was  large  enough  to  ensure  the  central  anomaly  would 
have  the  properties  of  a  global  solution,  we  simulated  the  recovery  of  a 
larger  15°x25°  region.  Our  results  actually  improved  slightly  for  the 
central  region,  giving  us  confidence  that  our  conclusions  would  be 
conservative. 

Another  problem  in  simulating  the  recovery  of  l°xl°  gravity  anomalies  is 
the  tremendous  amount  of  required  orbital  computations.  In  our  example 
using  a  region  15°  wide,  many  months  of  orbiting  by  polar  satellites  would 
have  been  required  to  obtain  dense  overflights  of  the  area.  In  addition,  a 
5-second  numerical  integration  step  size,  with  a  multistep  method  is  re- 
quired because  of  the  sharp,  high  frequency  effect  of  l°xl°  anomalies.  To 
avoid  this  lengthy  computation  we  considered  independent  orbital  arcs  only 
over  the  region  being  considered.  However,  the  initial  conditions  for  each 
arc  were  not  adjusted.  Our  assumption  was  that  the  ephemerides  of  the  two 
satellites  would  be  computed  in  long  arc  (a  few  days)  solutions  from  ground 
tracking. 

The  separation  of  the  ephemeris  computation  from  gravity  anomaly 
estimation  requires  justification,  especially  in  relation  to  the  effect  of 
systematic  errors  on  the  ephemerides.  For  example,  the  locations  of  ground 
tracking  stations  are  uncertain  to  some  degree,  and  this  error  propagates 
into  the  intersatellite  range-rate  measurement.  Fortunately,  this  effect  is 
readily  evaluated. 

The  usual  least-squares  solution  relates  observation  residuals,  80,  to 
adjusted  parameters,  a,  in  the  presence  of  data  noise,  e,  by 


10 


60  =  A6a  +  e  (1) 

where  A  is  a  matrix  of  partial  derivatives  of  the  measurements  with  respect 
to  the  adjusted  parameters.  If  unadjusted  parameters  with  error,  y,  related 
to  the  observations  by  a  matrix  of  partials,  K,  are  present,  then  the  obser- 
vation equation  becomes 

60  =  A  6a  -  K  y  +  e.  (2) 

If  W -1  is  the  covariance  matrix  of  the  measurement  noise,  e,  the  minimum 
variance  solution  6a,  which  ignores  the  error  in  the  unadjusted  parameters, 
is  biased  by  an  amount 

-(AT  W  A)"1  AT  W  K  y,  (3) 

The  effect  of  a  unit  value  of  the  y  parameters  is  given  by 

a(6a)/9y  =  -(AT  W  A)"1  AT  W  K.  (4) 

By  using  eq.  (4),  for  example,  we  can  predict  in  advance  the  effects  of  an 
error  in  station  location  on  any  particular  solution.  This  is  important  to 
evaluate,  especially  in  the  present  case,  because  of  the  smallness  of  the 
perturbation  of  a  l°xl°  mean  gravity  anomaly. 

Although  eq.  (4)  is  simple  in  principle,  it  is  computationally  complex. 
For  example,  we  can  easily  imagine  a  mission  where  10  Doppler  or  laser 
ground  tracking  stations  would  each  observe  two  satellites  for  6  months. 
The  partial  derivative  matrices  A  and  K  must  be  obtained  by  numerical  inte- 
gration of  variational  equations,  and  there  are  scores  of  these  equations. 

Many  sources  of  potential  systematic  error  can  arise  in  orbit  determina- 
tion. The  most  prominent  are  geopotential  and  tracking  station  coordinate 
errors.  We  did  not  consider  geopotential  error  here  because  we  have 
assumed,  of  course,  that  the  geopotential  was  being  estimated.  However, 
tracking-station  coordinate  errors  need  to  be  considered. 
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Using  the  computer  program  ORAN,  which  was  developed  around  eq.  (4),  we 
determined  the  effect  of  a  10-cm  uncertainty  in  the  height  of  one  station 
based  on  the  initial  conditions  and  the  intersatellite  range-rate  measure- 
ment, p,  for  two  satellites  spaced  3°  apart  in  the  same  200 -km  altitude 
polar  orbit.  The  orbital  arc  length  was  1  day  with  15  Doppler  tracking 
sites  well  -  distributed  around  the  globe,  using  line-of-sight  velocity 
measurements  of  a  =  1  cm/s  quality.  The  10-cm  value  for  the  error  of 
station  coordinates  is  realistic  for  the  1985  time  frame.  Figure  4  shows 
the  results. 

Two  significant  aspects  are  apparent  in  figure  4.  First,  notice  the  long 
period  of  the  oscillation  produced  by  the  station  coordinate  error.  This 
effect,  as  noted  before,  occurs  for  orbital  arcs  more  than  a  few  revolutions 
long  and  is  caused  by  the  dynamics  of  orbital  motion.  With  modern  instru- 
ments, the  timing  error  of  observations  is  negligible,  so  that  the  orbital 
period  is  extremely  well  -determined  in  long-arc  solutions  by  repeated 
passages  over  the  tracking  stations.  Thus  the  satellite  ephemeris  can  be  in 
error  only  periodically,  as  seen  in  figure  4.  The  effect  of  an  error  in  the 
location  of  an  observing  station  is  different  from  the  effect  of  a  localized 
gravity  anomaly,  such  as  a  l°xl°  block,  because  the  gravity  anomaly  causes  a 
distinct  local  change  in  the  acceleration  of  the  satellite. 

The  other  important  feature  of  figure  4  is  the  smallness  of  the  effect  of 
the  station  coordinate  error  on  the  intersatellite  range-rate.  The  ampli- 
tude is  2.5  x  10~4  cm/s,  about  5  percent  of  the  impulse-like  effect  of  a  10- 
mgal,  l°xl°  anomaly  block.  Because  the  wavelength  of  the  station-error 
induced  signal  is  more  than  10  times  the  wavelength  of  the  anomaly  block 
effect,  one  would  anticipate  little  aliasing  of  one  signal  into  the  other. 
But  it  is  a  simple  matter  to  filter  the  intersatellite  range-rate  signal  and 
eliminate  the  error  signal.  To  illustrate  this,  the  low  amplitude  curve  in 
figure  4  shows  the  effect  of  the  station  height  error  on  the  intersatellite 
range-rate  difference,  Ap  taken  over  1  minute.  Figure  5  shows  the  corre- 
sponding differences  of  the  range-rate  oscillation  produced  by  the  10-mgal, 
l°xl°  block.  As  expected,  the  long  wavelength  signal  is  greatly  reduced, 
and  the  high  frequency   impulse-like  signal  is  amplified.   In  terms  of 
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range-rate  differences,  the  ratio  of  the  10-mgal,  l°xl°  block  signal  to  a 
station  coordinate  error  signal  is  about  400.  The  conclusion  reached  from 
this  discussion  is  that  we  do  not  need  to  consider  further  the  effect  of 
error  in  tracking  station  coordinates  on  our  simulation  of  the  recovery  of 
gravity  anomalies.  Such  errors,  including  tracking-data  systematic  errors 
that  also  produce  long  wavelength  ephemeris  errors,  are  readily  filtered 
from  the  basic  data  type,  the  intersatellite  range-rate.  Our  simulations, 
therefore,  consider  only  the  contribution  of  the  statistical  uncertainty  of 
the  intersatellite  range-rate  data  to  the  uncertainty  of  the  estimated  l°xl° 
mean  gravity  anomalies.  We  filtered  the  signal  by  taking  successive  differ- 
ences. Incidentally,  this  simple  filter  was  merely  a  convenience.  Much 
better  filters  can  be  constructed. 

SIMULATION  RESULTS 

We  have  seen  that  the  perturbation  of  the  intersatellite  velocity  depends 
on  the  relative  separation  of  the  satellites.  Figure  3  shows  that  the 
perturbation  is  maximized  for  l°xl°  blocks  and  200-km  altitude  by  a  separa- 
tion exceeding  5°.  However,  the  accuracy  of  the  intersatellite  tracking 
device  itself  depends  on  the  separation  of  the  satellites.  We  assumed 
intersatellite  tracking  of  10 "^  mm/s  precision  with  an  integration  time 
(i.e.,  data  interval)  of  8  seconds.  According  to  estimates  made  at  The 
Johns  Hopkins  University/Applied  Physics  Laboratory,  this  rate  and  precision 
are  obtainable  with  a  microwave  device  for  a  separation  of  3°  or  less 
(Fischell  1979).  Figure  3  shows  that  this  separation  does  not  give  the 
maximum  amplitude  intersatellite  range-rate  signal.  The  maximum  possible 
range-rate  signal  is  about  50  percent  greater. 

The  region  chosen  for  our  initial  simulation  experiments  was  15°xl5°, 
containing  225  unknown  mean  anomalies,  and  centered  at  0°  latitude  and 
longitude.  Figure  6  shows  the  upper  right  side  (NE  quadrant)  of  this  region 
and  the  geometry  of  the  downward  satellite  tracks.  Six  downward  tracks  were 
selected  over  each  1°  for  the  central  4°  of  longitude  and  three  tracks  over 
each  subsequent  degree  giving  a  total  of  64  downward  tracks  over  the  region. 
Figure  6  shows  that  the  Earth's  rotation  has  the  effect  of  causing  the 
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ground  tracks  of  the  polar  orbiting  satellite  to  be  slightly  inclined  to  a 
meridian.  Note  that  the  right  edge  of  the  region  is  carried  away  from  the 
tracks,  resulting  in  poor  sampling.  This  is  one  of  the  reasons  our  results 
are  valid  only  for  the  central  area  of  the  15°xl5°  region. 

Our  first  simulation  (fig.  7)  shows  the  results  for  the  case  of  64  arcs  of 
the  two  polar  orbiting  satellites  at  200-km  altitude,  spaced  3°  apart,  and 
passing  in  the  north-south  direction  over  the  region.  The  numbers  in  each 
square  (many  are  omitted  for  clarity)  are  the  standard  deviations  of  the  mean 
gravity  anomalies  in  milligals  for  each  block.  These  were  taken  from  the 
computed  covariance  matrix  of  the  solution,  assuming  an  intersatellite 
range-rate  measurement  precision  of  10 "®  mm/s,  but  filtered  by  taking  1  min- 
ute differences  every  8  seconds.  The  results  for  the  central  region  are 
about  6  mgal. 

Figure  8  shows  the  standard  deviations  in  milligals  for  the  l°xl°  mean 
anomalies  after  adding  to  the  simulation  in  figure  7  the  equivalent  64 
north-going  passes  of  the  two  satellites.  This  128-arc  simulation  is  im- 
proved by  a  factor  of  2  rather  than  the  40  percent  improvement  expected  by 
doubling  the  amount  of  data. 

The  128-arc  simulation  gives  12  passes  over  each  of  the  central  6°  of 
longitude  with  the  remaining  passes  distributed  evenly  over  the  outer 
blocks.  This  reduction  for  the  outer  blocks  was  done  only  for  economy  and, 
if  anything,  makes  the  results  conservative. 

The  128-arc  simulation  is  encouraging  because  the  standard  deviations  for 
the  central  region  are  about  3  mgal.  The  peripheral-region  results  are 
better,  but  only  because  the  blocks  outside  the  15°xl5°  region  are  not  ad- 
justed. Again,  the  peripheral  block  a's  are  not  realistic.  In  addition, 
the  right-side  edge  values  are  larger  than  the  left-side  values.  As  noted 
previously,  this  occurs  because  the  rotation  of  the  Earth  carries  the  region 
away  from  the  satellite  passes. 

Of  equal  importance  to  the  a's  are  the  correlation  coefficients.  We 
anticipated  high  values  because  the  perturbation  caused  by  a  l°xl°  block 
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endures  for  a  relatively  long  scale  compared  to  the  block  size.  Figure  9 
shows  the  correlation  coefficients  of  the  center  block  with  the  other  blocks 
for  the  simulation  of  figure  8.  The  correlation  of  the  center  block  with 
its  east  or  west  neighbor  is  very  high  (-0.88).  Thus,  the  block  estimation 
errors  caused  by  data  noise  will  have  almost  the  same  magnitude  of  error  but 
in  an  opposite  direction  to  its  immediate  neighbor  in  the  east-west  direc- 
tion. We  should  expect,  therefore,  that  the  results  for  2°x2°,  or  larger, 
blocks  would  be  better.   This  is  confirmed  below. 

Considering  the  still  significant  correlations  at  the  north-south  borders 
in  figure  9  compared  to  the  smaller  values  at  the  east-west  borders,  we  need 
to  consider  now  whether  or  not  the  overall  region  of  15°xl5°  was  really 
large  enough  for  the  center  blocks  to  have  the  statistical  properties  of  a 
global  solution.  Figure  10  shows  the  results  for  a  64-arc  (north-to-south 
only)  simulation  of  the  solution  for  375  l°xl°  blocks  with  the  orbital  arcs 
and  the  region  extended  from  15°xl5°  to  15°x25°.  (This  limited  case  of  one 
direction  was  chosen  for  economy.)  We  can  compare  this  directly  with  the 
equivalent  64-arc  case  of  225  l°xl°  blocks  shown  in  figure  6.  We  see  that 
the  15°x25°  case  is  actually  improved  by  one-third  as  compared  to  the 
15°xl5°  situation.  Recalling  that  the  128-arc  15°xl5°  combined  solution, 
which  used  both  the  north-going  and  south-going  passes,  was  twice  as  good  as 
either  the  north-  or  south-going  solution  independently,  we  can  reasonably 
conclude  that  a  complete  128-arc,  15°x25°  region  solution  would  be  about 
2  mgal  for  the  uncertainty  of  the  blocks  in  the  central  area.  However,  the 
east-west  correlations  are  not  significantly  improved. 

Having  obtained  a  basic  result  of  a  few  milligals  for  the  200-km  altitude,  we 
now  consider  the  impact  for  missions  at  other  altitudes.  Figure  11  shows 
the  a's  for  the  128-arc  case  (considered  previously  in  fig.  8)  with  the 
orbiting  pair  instead  at  150-km  altitude.  Figure  12  gives  the  results  for 
an  altitude  of  250  km.  The  150-km  case  is  better  by  about  a  factor  of  5  and 
the  250-km  case  is  worse  by  about  a  factor  of  5  than  the  200-km  situation. 
In  addition,  the  correlation  with  the  neighboring  block  drops  to  0.81  for 
the  150-km  example  and  rises  to  0.92  for  the  250-km  altitude.  The  sharp 
improvement  obtained  by  decreasing  the  altitude  toward  the  scale  of  the 
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l°xl°  blocks  (-111  km/deg  at  the  surface)  tends  to  confirm  Schwarz's  (1971) 
rule  that  the  horizontal  resolution  of  the  low-low  satellite  experiment  is 
approximately  equal  to  the  altitude  of  the  orbiting  pair. 

Our  final  consideration  was  the  possible  performance  of  a  so-called 
"high-low"  system  where  one  of  the  satellites  is  geosynchronous.  Several 
advantages  of  this  method  have  been  alleged.  Indeed,  for  a  given  intersatel- 
lite  speed  measurement  specification,  the  high-low  case  should  be  superior 
overall  because,  for  blocks  away  from  the  Equator,  both  the  horizontal  and 
vertical  components  of  the  intersatellite  velocity  perturbation  are  observed. 
Figure  13  gives  the  results  obtained  for  a  64-arc  case  (south-going  passes 
alone)  with  the  geosynchronous  member  of  the  pair  placed  directly  over  the 
center  of  the  15°xl5°  region  and  the  low  satellite  at  150-km  altitude.  Com- 
parison with  f  igure  10  shows  that  this  high-low  configuration  of  64-arcs  is 
as  effective  as  a  128-arc  low-low  experiment  with  all  else  being  equal.  But 
we  emphasize  these  results  depend  entirely  on  intersatellite  range-rate 
measurements  of  10  3  mm/s  precision.  This  precision  appears  to  be  obtainable 
with  laser  interferometric  or  sophisticated  microwave  devices  for  satellites 
orbiting  close  to  each  other,  but  probably  further  technological  development 
is  required  for  observations  to  ,  a  geosynchronous  satellite  instead  of  an- 
other near-Earth  one. 


It  is  interesting  to  compute  the  results  obtainable  when  the  block  size  is 
changed  to  2°x2°.  If  the  l°xl°  block  errors  were  uncorrelated,  mean  2°x2° 
anomalies  would  have  a  standard  deviation  of  \  of  the  l°xl°  anomalies 
comprising  them.  However,  because  of  the  large  negative  correlation  of  ad- 
jacent blocks,  2°x2°  mean  anomalies  have  a  precision  of  0.29  mgal,  an  im- 
provement of  10  compared  to  the  l°xl°  case.  Because  the  200-km  altitude 
corresponds  almost  exactly  to  the  dimension  of  a  2°x2°  block,  this  result 
provides  another  demonstration  of  Schwarz's  rule. 

DISCUSSION 

It  would  be  desirable  to  be  able  to  give  a  firm  number  for  the  capability 
of  an  intersatellite  tracking  mission  to  determine  the  geopotential,  but  our 
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study  did  not  achieve  this.  However,  from  a  conservative  viewpoint  we  be- 
lieve that  we  have  answered  an  important  question:  Is  the  capability  of 
satellite-to-satellite  tracking  for  determining  the  geopotential  at  the 
l°xl°-resolution  level  near  10  mgal  or  1  mgal?  Our  answer  is  the  latter, 
especially  because  our  simulations  were  so  conservative.  For  example,  the 
number  of  passes  over  the  blocks  in  our  simulation  would  be  attained  in  a 
real  mission  in  about  4  months.  Much  longer  missions  are  within  the  capa- 
bility of  surface-force  compensation  systems,  so  better  results  could  be 
expected  from  longer  missions.  In  addition,  a  somewhat  lower  altitude  could 
be  chosen.  Recall  that  dropping  from  200  km  to  150  km  caused  a  five-fold 
improvement.  An  intermediate  altitude  choice  of  about  180  km,  rather  than 
200  km,  could  result  in  a  substantial  improvement,  but  at  the  cost  of  a  50 
percent  increase  in  atmospheric  drag  force.  Improvements  are  also  possible 
in  intersatellite  measurement  capabilities.  This  would  enable  an  increase 
in  the  intersatellite  separation  and  would  improve  the  signal/  noise  for  the 
experiment.  Finally,  a  careful  inclusion  of  a  priori  information  taken  from 
the  vast  store  of  existing  geopotential  data  would  improve  the  situation. 
Thus,  it  seems  apparent  that  a  suitably  optimized  1  year  low-low  satellite 
mission  could  produce  mean  anomalies  at  the  1-mgal  level  of  precision.  The 
high  correlations  in  the  east- west  direction  are  unfortunately  large,  but 
are  tolerable  because  the  a's  are  so  small  compared  to  the  actual  variation 
of  gravity  anomalies  (about  30  mgal  rms  in  l°xl°  blocks) . 

Some  remarks  should  be  made  on  the  significance  of  a  global  l°xl°  precise 
gravity  field.  Such  a  field  does  not  necessarily  produce  an  ocean  geoid 
suitable  as  a  reference  surface  for  satellite  altimetry  at  the  1°  wavelength. 
Chovitz  (1973)  showed  that  Kaula's  (1966)  rule-of-thumb  10~5/£2  for  the 
degree  variances  of  geopotential  coefficients  implies  unaccounted-for  geoid 
undulations  of  64/£  meters  for  a  geoid  computed  to  degree  £.  For  a  l°xl° 
solution,  the  undulations  shorter  than  1°  could  be  about  30  cm.  At  wave- 
lengths longer  than  1°,  the  geoid  obtained  from  l°xl°  mean  anomalies  would 
be  highly  accurate.  In  addition,  the  existence  of  a  geopotential  field 
accurate  up  to  the  l°xl°  scale  would  effectively  solve  the  problem  of  satel- 
lite orbit  determination.   With  such  a  field  and  a  satellite  equipped  with  a 
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surface-force  compensation  device,  a  satellite  ephemeris  accurate  to  better 
than  10  cm  could  be  obtained.  Also,  the  computation  of  gravimetric  geoids 
such  as  the  one  by  Marsh  and  Chang  (1978),  now  widely  used  as  a  reference 
surface  for  satellite  altimetry,  could  be  improved  because  the  effects  of 
mean  anomalies  in  regions  distant  from  a  geoid  computation  point  would  be 
accurately  known.  Finally,  a  global  l°xl°  field  would  be  of  interest  to 
geophysicists . 
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